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Abstract. We present an approach for analyzing initial-boundary value problems 
for integrable equations whose Lax pairs involve 3x3 matrices. Whereas initial 
value problems for integrable equations can be analyzed by means of the classical 
Inverse Scattering Transform (1ST), the presence of a boundary presents new chal- 
lenges. Over the last fifteen years, an extension of the 1ST formalism developed by 
, Fokas and his collaborators has been successful in analyzing boundary value prob- 

(""*) lems for several of the most important integrable equations with 2x2 Lax pairs, 

CNj such as the Korteweg-de Vries, the nonlinear Schrodinger, and the sine-Gordon 

^ equations. In this paper, we extend these ideas to the case of equations with Lax 

pairs involving 3x3 matrices. 
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1. Introduction 

Several of the most important PDEs in mathematics and physics are integrable. 
Integrable PDEs can be analyzed by means of the Inverse Scattering Transform (1ST) 
formalism, whose introduction was one of the most important developments in the 
theory of nonlinear PDEs in the 20th century. Until the 1990's the 1ST methodol- 
ogy was pursued almost entirely for pure initial value problems. However, in many 
laboratory and field situations, the wave motion is initiated by what corresponds to 
£N) the imposition of boundary conditions rather than initial conditions. This naturally 

qq leads to the formulation of an initial-boundary value (IBV) problem instead of a pure 

initial value problem. 

In 1997, Fokas announced a new unified approach for the analysis of IBV problems 
for linear and nonlinear integrable PDEs [12|, LT3] (see also |14j ) . The Fokas method 
provides a generalization of the 1ST formalism from initial value to IBV problems, 
^ and over the last fifteen years, this method has been used to analyze boundary value 

problems for several of the most important integrable equations with 2x2 Lax pairs, 
such as the Korteweg-de Vries, the nonlinear Schrodinger, the sine-Gordon, and the 
stationary axisymmetric Einstein equations, see e.g. [21 El El H31 HH [TBI EE EU [22], 
ESI ES]. Just like the 1ST on the line, the unified method yields an expression for 
the solution u(x, t) of an IBV problem in terms of the solution of a Riemann-Hilbert 
(RH) problem. In particular, the asymptotic behavior of u(x, t) can be analyzed in 
an effective way by using this RH problem and by employing the nonlinear version of 
the steepest descent method introduced by Deift and Zhou [1 1J . 

The purpose of this paper is to develop a methodology for analyzing IBV problems 
for integrable evolution equations with Lax pairs involving 3x3 matrices. Although 
the transition from 2 x 2 to 3 x 3 matrix Lax pairs involves a number of novelties, the 
two main steps of the method of |121 113] remain the same: (a) Construct an integral 
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representation of the solution characterized via a matrix RH problem formulated 
in the complex /c-plane, where k denotes the spectral parameter of the Lax pair. 
This representation involves, in general, some unknown boundary values, thus the 
solution formula is not yet effective, (b) Characterize the unknown boundary values 
by analyzing the so-called global relation. In general, the characterization of the 
unknown boundary values involves the solution of a nonlinear problem. However, for 
so-called linearizable boundary conditions, this problem can be by-passed since the 
unknown boundary values can be eliminated using only algebraic manipulations. 

In this paper, we will show how steps (a) and (b) can be implemented for a proto- 
typical example of an equation with a 3 x 3 Lax pair. We expect that a similar analysis 
will apply also to other integrable equations with 3x3 Lax pairs, both on the half-line 
and on the interval, although in certain cases additional technical difficulties will arise. 
For example, for the Boussinesq equation the behavior of the eigenfunctions near the 
point k = requires special attention^] Physically relevant equations with 3x3 Lax 
pairs include the Boussinesq, Degasperis-Procesi [9], Kaup-Kupershmidt [20], Sasa- 
Satsuma |27j . Sawada-Kotera |28j . two-component vector nonlinear Schrodinger [23], 
and 3-wave resonant interaction [31] equations. 

1.1. The transition from 2 x 2 to 3 x 3 Lax pairs. Let us comment on some of 
the implications of the transition from 2 x 2 to 3 x 3 Lax pairs. The implementation 
of the step (a) mentioned earlier in the case of 2 x 2 matrix Lax pairs is achieved 
by introducing eigenfunctions {fij(x, t, k)} which are defined by integration from the 
'corners' of the physical domain, see [13] • The column vectors of the /ij's are bounded 
and analytic in different sectors of the complex fc-plane and these column vectors 
are easily combined into a sectionally analytic function suitable for the formulation 
of a RH problem. However, for a Lax pair involving 3 x 3-matrices, the bounded 
and analytic eigenfunctions suitable for the formulation of a RH problem involve 
rather complicated combinations of the entries of the /Uj's. Moreover, because of 
the limited domains of boundedness of the /i/s, the boundedness properties of these 
combinations are not evident. We will therefore use a different approach for finding 
these combinations and their boundedness domains. Instead of taking the fij's as our 
starting point, we will define analytic eigenfunctions, denoted by {M n (x, t, k)}, via 
integral equations which involve integration from all three corners simultaneously. In 
the absence of bound states, this formulation is adequate. However, since the integral 
equations defining the eigenfunctions now are of Fredholm rather than Volterra type, 
there may exist points {kj}, kj 6 C, at which the M n 's have singularities. In order to 
deal with these singularities (which are related to the existence of solitons), we will 
relate the M n 's to the //j's by solving a matrix factorization problem. 

Another difference in the implementation of step (a) is that the RH problem in the 
case of a 2 x 2 Lax pair splits the complex fc-plane into 4 sectors, whereas a larger 
number of sectors is in general required in the case of 3 x 3 Lax pairs, e.g. for our 
main example, the RH problem splits the complex fc-plane into 12 sectors. 

The implementation of the step (b) mentioned earlier involves eliminating the 
unknown boundary values from the formulation of the RH problem. In the case of 
linearizable boundary conditions for equations with 2x2 matrices, this elimination is 
achieved by algebraic manipulation of the so-called global relation and the equations 
obtained from the global relation under certain transformations in the /c-plane. As 
shown in section [4] below, similar ideas can be used to analyze linearizable boundary 
conditions in the case of 3 x 3 Lax pairs. However, in the 3x3 case the algebraic 
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relations cannot always be directly used to eliminate the unknown boundary values 
from the definition of the jump matrices. Instead, we will first analytically extend 
the domain of definition of the jump matrix, before we utilize the algebraic relations. 
Finally, we perform another analytic continuation to find the expression for the jump 
matrix on the relevant contour. The analyticity domains of the involved matrices are 
just sufficient to allow for this approach. 

Let us finally point out that some pioneering works on the theory of inverse scatter- 
ing for initial- value problems for equations with 3x3 Lax pairs are [U [TO] 120] . Other 
examples of the use of piecewise analytic solutions in studying the inverse problems 
for n x n systems (particularly for n = 4) can be found in (3J HI [29]. 

1.2. Organization of the paper. Our main example is introduced in subsection 
1.3 The spectral analysis of the associated Lax pair is performed in section [2} In 



section [3j we formulate the main RH problem, and this concludes the implementation 
of step (a) above. 

Step (b) is implemented in sections [4] and [5] In section |4j we consider linearizable 
boundary conditions, whereas nonlinearizable boundary conditions are analyzed in 
section [5} Section [6] contains some concluding remarks. In appendix [Aj we explain 
the relationship between the constructions of this paper and the formalism of [121 IT3] 
for 2x2 Lax pairs. In appendix [Bl we use an extension of the standard Fredholm 
theory to study a set of integral equations. 

1.3. The main example. To be concrete, we will consider the system 



Ht + 75 + 2irr x = 0, 
ir t - \r xx + 2iqq x = 0, 



where q(x,t) and r{x,t) are complex- valued functions of (x,t) £ f2, with f2 denoting 
the half-line domain 

n = {0 < x < oo, < t < T} 

and T > being a fixed final time. This system is the compatibility condition of the 
3x3 Lax pair 

'\ N 

J = I u? | , uj = eT, (1.2) 



ip x - kJijj = Viip, r _f n (i2 J ,. _?p 

f/H + k 2 J 2 ^ = V 2 ip, 





where ^(x,t,k) is a 3 x 3-matrix valued function, k £ C is the spectral parameter, 
and {Vi(x, t), V 2 {x, i, k)} are 3 x 3-matrix valued functions given by 

V 2 (x,t,k) = kV^ ] (x,t) + vj 0) (x,t), (1.3) 

UJ 2 r\ . ( q x -r x \ / r 2 q 2S 

Q , W 0) = 4 ~r x q x - \q 2 r 2 
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We will denote the initial data of (1.1) by {qo(x), ro(x)}, while the Dirichlet and 
Neumann boundary values will be denoted by {go(t), ho(t)} and {gi(t), hi(t)}, re- 
spectively, i.e. 

q(x,0) = qo(x), r(x,0) = ro(x), < x < oo; 

q(0,t) = g (t), r(0,t) = /*,(*), (1.4) 

Qx(0,t) = gi (t), r x {0,t) = hx(t), 0<t<T. 
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Remark 1.1. 1. We have chosen (1.1) as our main example, because the Lax pair 
( 1.2 ) is a natural 3x3 generalization of the 2x2 Lax pair for the nonlinear Schrodinger 
equation. Indeed, the x and t parts of the Lax pair for NLS involve the matrices ka^ 
and k 2 <T3, respectively, where 03 is the diagonal matrix whose entries are the second 
roots of unity, i.e. 03 = diag(l, —1). Analogously, the x and t parts of (1.2) involve 
the matrices kJ and A; 2 J 2 , where J is the diagonal matrix whose entries are the third 
roots of unity. 

2. When r = q, the system (1.1) reduces to the following integrable equation (cf. 
Eq. (1.5) in pU): 

1 . . 

iq t + -j=q xx + 2iqq x = 0. (1.5) 

2. Spectral analysis 

Our goal in this section is to define sectionally analytic eigenfunctions of the Lax 
pair (1.2) which are suitable for the formulation of a Riemann-Hilbert problem. 

2.1. The closed one-form. Suppose that q(x,t) and r(x,t) are sufficiently smooth 
functions of (x, t) in the half-line domain f2 which decay as x — > 00. Introducing 

L(x,t,k) and Z(x,t,k) by 

L = kJ + Vi, Z = -k 2 J 2 + V 2 , (2.1) 



(2.2) 



we can write (1.2) as 

\i(j x = Lip, 

(> = Zip. 

The functions L and Z satisfy 

trL = 0, trZ = 0, 

as well as the Z3 symmetry 



L(k) = AL(uk)A-\ Z(k) = AZ(uk)A'\ A=\0 1 . (2.3) 




In order to keep the notation concise, we define diagonal matrices C and Z by 
C = lim L = kJ, Z = lim Z = -k 2 J 2 , 



X—I-OO 

|3 „„J/,./i-n3. 



and denote the diagonal entries of these matrices by {lj{k)}f &nd{zj(k)}f : 

C = diag(/i,/ 2 ,^), Z = diag(zi, 212,213). 
Introducing a new eigenfunction /i(x, t, k) by 

^ = ne Cx+zt , (2.4) 



we find the Lax pair equations 



fj, x - [C,/Jt] = Vifj,, 



(2.5) 



Letting C and Z denote the operato rs wh ich act on a 3 x 3 matrix X by CX = [C, X] 
and ZX = [Z,X], the equations in (2.5) can be written in differential form as 



d 



(e- ia; -*V) = W, (2-6) 
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Figure 1. The contours 71, 72, and 73 m i/ie (x,t)-plane. 



73 



where W(x, i, fc) is the closed one- form defined by 

t (V 1 dx + V 2 dt) f x. 



W = e - Cx - zt i 



Remark 2.1. In the case of equation (1.5), the assumption q 



Lax pair possesses the following additional ^-symmetry: 



(2.7) 

= r implies that the 

A o\ 

L{k) = BL(k)B , Z{k) = BZ{k)B~ 1 , B= 1 . 

V° 1 °/ 

2.2. The M n 's. Let {7j}f be contours in the (x,t)-plane which connect (xj,tj) to 
(x,t), where 

(xi, ti) = (0, T), (x 2 , t 2 ) = (0, 0), (x 3 , i 3 ) = (00, t). 

We will assume that T < 00, except for in section [4] where it will be assumed that 
T = 00. Choosing the particular contours shown in Figure [TJ we have the following 
inequalities on the contours: 



71 : x - x' > 0, i - t' < 0, 

72 : x — x > 0, i — i' > 0, 

73 : x — x < 0. 



(2. 



For each n = 1, . . . , 12, define a solution M n (x, t, k) of (2.5) by the following system 
of integral equations: 

(M n ) ij (x,t,k) = 6 ij + I (e^+^Wnix'J,^)^, keD n , i,j = 1,2,3, 

(2.9) 

where W n is given by (|2.7|) with \i replaced by M n , the contours 7^, n = 1, . . . , 12, 



i,j = 1,2,3, are defined by 

71 if Keli(k) < Kelj(k) and Rezi(k) > Kezj(k) 
iij = { 72 if Rek(k) < Relj(k) and Re Zi(k) < Re zj (k) for k e D n , 
73 if ReZi(fc) > ReZj-(fc) 

(2.10) 
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and {.D n }J 2 denote twelve open, pairwisely disjoint subsets of the complex Zc-plane 
defined by (see Figure [2]) 

Di = {k G C I Re Z3 < Re I2 < Re l\ and Re z\ < Re z% < Re Z2}, 

D2 = {k G C I Re Z3 < Re I2 < Re l\ and Re z% < Re z\ < Re Z2}, 

-D3 = {k G C I Re Z3 < Re l\ < Re I2 and Re z% < Re Z2 < Re zi}, 

D4 = {k G C I Re Z3 < Re l\ < Re I2 and Re Z2 < Re zj, < Re z{\, 

D5 = {/c G C I Re Zi < Re Z3 < Re I2 and Re Z2 < Re z\ < Re 23}, 

L>6 = £ C I Re Zi < Re Z3 < Re h and Re z\ < Re 2:2 < Re 2:3}, 

Dj = {k G C I Re l\ < Re Z2 < Re Z3 and Re z\ < Re 23 < Re Z2}, 

D$ = {k G C I Re Zi < Re I2 < Re Z3 and Re z% < Re z\ < Re Z2}, 

Dq = {k G C I Re I2 < Re l\ < Re Z3 and Re z% < Re Z2 < Re zi}, 

D10 = {k G C I Re I2 < Re l\ < Re Z3 and Re ^2 < Re Z3 < Re zi}, 

D\\ = {k G C I Re Z2 < Re Z3 < Re l\ and Re Z2 < Re zi < Re 2:3}, 

D12 = {k G C I Re Z2 < Re Z3 < Re l\ and Re z\ < Re Z2 < Re Z3}. 

The following proposition ascertains that the M n 's defined in this way have the prop- 
erties required for the formulation of a RH problem. 

Proposition 2.2. For each n = 1, . . . , 12, the function M n (x, t, k) is well-defined by 



equation (2.9) for k G D n and (x,t) G For any fixed point (x,t), M n is bounded 
and analytic as a function of k G D n away from a possible discrete set of singularities 
{kj} at which the Fredholm determinant vanishes. Moreover, M n admits a bounded 
and continuous extension to D n and 

M n (x,t,k) = I + 0(l/k), k^oo, k£D n . (2.11) 
Proof. The boundedness and analyticity properties are established in appendix lB| 



here we prove (2.11). Substituting the expansion 



M = + — + -p- + • • • , k -> 00, 



into the x-part of the Lax pair (2.5), the terms of O(k) yield 



[J,M (0) ] = i.e. M (0) is a diagonal matrix. 



The diagonal terms of 0(1) imply that M^ ' is a constant matrix. Evaluation of (2.9) 



at (x 3 ,t 3 ) implies that = I. □ 

2.3. The jump matrices. We define spectral functions S n (k), n = 1, . . . , 12, by 

S n (k) = M n (0,0,k), k£D n , n = l,...,12. (2.12) 

Let M denote the sectionally analytic function on the Riemann Zc-sphere which equals 
M n for k G D n . Then M satisfies the jump conditions 

M n = M m J m ^ n , k£D n C]D m , n,m = 1, ...,12, n ^ m, (2.13) 

where the jump matrices J m ,n{x, t, k) are defined by 

Jm,n = e tx+± \S-}S n ), n,mG{l,...,12}. (2.14) 



Since the integral equations (2.9) that define M n (0, 0, k) only involve integration along 



the initial half-line {0 < x < 00, t = 0} and along the boundary {x = 0, < t < T}, 
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Figure 2. The sets D n , n = 12, which decompose the complex k-plane. 



the SVi's (and hence also the J m ,n's) can be computed from the initial and bound- 



ary data alone. Thus, relation ( 2.13| ) provides the jump condition for a RH prob- 
lem, which, in the absence of singularities, can be used to reconstruct the solution 
{q(x, t), r(x, t)} from the initial and boundary data. However, if the M n 's have pole 
singularities at some points {kj}, kj £ C, the RH problem needs to include the 
residue conditions at these points. For the purpose of determining the correct residue 
conditions (and also for the purposes of analyzing the linearizable and nonlineariz- 
able boundary conditions in sections [4] and [5]) , it is convenient to introduce three 
eigenfunctions {fij(x,t, k)}f in addition to the M n 's. 



2.4. The /%-'s. We define three eigenfunctions {/%}? of (2.5) by the Volterra integral 
equations 

Hj{x,t,k)=I+ f e t{ - k)x+ ^ k)t Wj{x',t',k), j = 1,2,3, (2.15) 



where Wj is given by (2.7) with \x replaced with The third column of the matrix 
equation (2.15) involves the exponentials 



,(h-l 3 )(x-x')+(zi-z 3 )(t-t') 



,(h-l3)(x-x')+(z 2 -z 3 )(t-t') 



Using the inequalities in (2.8) it follows that these exponentials are bounded in the 



following regions of the complex /c-plane: 

{ReZi < Re/ 3 } n {Re/ 2 < Re/3} n {Rez 3 < Rezi} n {Rez 3 < Rez 2 }, 
{Reh < Rel 3 } n {Re/ 2 < Re/3} n {Rezi < Rez 3 } n {Rez 2 < Re2 3 }, 
{Re/ 3 < Reh} D {Re/ 3 < Re/ 2 }. 



71 

72 
73 



Since the equations in ( 2.15| ) are Volterra integral equations, these boundedness prop- 
erties imply that the third column vectors of [i\ and /x 3 are bounded and analytic 
for k G C such that k belongs to 1Z and S, respectively, where the sets 1Z and S are 
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Figure 3. The decompositions of the complex k-plane defined by the sets 7Z, S, T 
and their reflections 1Z, S, T . 

defined by (see Figure [3J 

K = D 8 \JD 9 , 5 = DiUD 2 UD 3 UD 4 . 
Similar conditions are valid for the other column vectors. Thus, 

t, k) is bounded for k G (ujTZ, uj 2 1Z, 7Z), 
mix, t, k) is bounded for k G (ujS, uj 2 S, S), 

where k £ (uj7Z,ijj 2 1Z,7Z) indicates that the first, second, and third columns of the 
equation are valid for k in the sets ujTZ, lo 2 1Z, and 1Z, respectively. On the other hand, 
the corresponding sets where the columns of (12 are bounded are all empty. 

Note that /Ui and fi2 are entire functions of k, whereas the first, second, and third 
columns of [13, whose definitions involve integration from x = 00, are defined only for 
k G (oaS, uj 2 S, S). Moreover, as k — > 00, 

m(x,t,k) = I + 0(l/k), k^oo, k G (ujU,uj 2 TZ,TZ), 
fj, 3 (x,t,k) = l + 0(l/k), k^oo, k G (loS, uj 2 S, S). 
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For x = 0, \i\ and [i 2 have the following enlarged domains of boundedness: 
/xi (0, t, k) is bounded for k G (ujTZ U ujTZ, lu 2 K U a; 2 7t, KUK), 
/x 2 (0, i, fc) is bounded for fc G (wT U ut, u?T U w 2 T, Tut), 

where T = D^UDq and, for a set U C C, U denotes the reflection of U with respect to 
the origin, i.e. U = {k G C | — k G U}; in particular, 1Z = D 2 UD 3 and T = -D11U.D12. 

Remark 2.3. We have defined two sets of eigenfunctions: {M n }} 2 and The 
unified approach of [12] for Lax pairs involving 2 x 2-matrices also implicitly utilizes 
two types of eigenfunctions: the fij's are used for the spectral analysis, whereas the 
RH problem is formulated in terms of another set of eigenfunctions; our M n 's are the 
analogs of this latter set of eigenfunctions, see appendix [A) 

2.5. The minors of the eigenfunctions. We will also need the analyticity and 
boundedness properties of the minors of the matrices fij(x,t,k), j = 1,2,3. To this 
end, we recall that the cofactor matrix X A of a 3 x 3 matrix X is defined by 

/ mn(X) -m l2 (X) m 13 (X) 
X A = -m 21 (X) m 22 (X) -m 23 (X) 
\ m 3 i(X) -m 32 (X) m 33 (X) 

where rriij(X) denotes the (ij)th minor X, i.e. m.ij(X) equals the determinant of the 
2 x 2-matrix obtai ned from X by deleting the ith row and the jth column. 
It follows from (2.5) that the eigenfunction jx satisfies the Lax pair 



A , [7 „Ai _ \rT ,,A ( 2 - 16 ) 



Thus, the eigenfunctions {/J-f}i are solutions of the integral equations 

fif(x,t,k) = I - f e-^-^-^-^iVfdx + V^dt)^, j = 1,2,3. (2.17) 

The functions fi A and \i 2 are entire functions of k, whereas the first, second, and 
third columns of /j a , whose definitions involve integration from x = 00, are defined 
for k G (u!<S,u 2 S,<S). Moreover, 

fi 2 (x,t,k) is bounded for k G (ujit,J 2 Tl,'lZ) 1 

jj, A (x,t,k) is bounded for k G (uaS, uj 2 S, S), 

whereas the corresponding sets where the columns of fi A are bounded are all empty. 
It also holds that 

fj, A (x,t,k) = l + 0(l/k), k^oo, k G (uTZ, u 2 1Z, TZ), 

l4(x,t,k) = I + 0(l/k), k^oo, k e (ooS,lj 2 S,S). 

For x = 0, fi A and fx 2 have the following enlarged domains of boundedness: 

fj,f(Q, t, k) is bounded for k G (wT U wT, to 2 T U w 2 T, T UT), 

fj,f(Q,t,k) is bounded for k G (uTZ U uTZ, u?1Z U oj 2 K, K U TZ). 

Remark 2.4. 1. The inverse of a unit determinant 3x3 matrix X is given by X^ 1 = 
(X A ) T . The tracelessness of the matrices L and Z implies that det Hj(x, t, k) = 1, 
k G C, for j = 1,2. 
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2. It is crucial that many of the minors m>ij(nk) have much larger domains of 
definitions and boundedness than the analysis in subsection 2.4 suggests. For exam- 
ple, we have seen that \ji3\2 and [/a^s are defined and bounded for k in uj 2 S and S, 
respectively, where [Xh denotes the jth column of the matrix X. Nevertheless, the 
combination mn(/x 3 ) = (^3)22(^3)33 - (^3)23(^3)32 is bounded not only for k > but 
for all k G lvS. 



2.6. Symmetries. The Z 3 symmetry (2.3) implies corresponding symmetries for the 
eigenfunctions. 

Lemma 2.5. The functions M and {^j}f obey the symmetries 
M(x,t,k) = AM(x,t,ujk)A~ l , k £ C, 

and 

fij(x,t,k) = Afij(x,t,u>k)A~ 1 , fceC if j = 1,2; 

k G (ojS,uj 2 S,S) if j = 3. 



(2.18) 



(2.19) 



Proof. We will prove (2.18) for k G D\\ the proofs when k G Dj for j / 1 are 
analogous. Thus, suppose that k G D±. Then tok G -D5 and the functions 

1>{x,t,k) = M!(x,t,A:)e £(fc)x,+2(fc) * and <j>{x,t,k) = M 5 (x,t,uk)e c(ujk)x+Z{uk)t 

satisfy the equations 



ipx(x, t, k) = L(x, t, k)if)(x, t, k), 
ipt(x, t, k) = Z(x, t, k)ijj(x, t, k), 



(j) x (x, t, k) = L(x, t, uk)(j)(x, t, k), 
4>t(x, t, k) = Z(x, t, uk)(j)(x, t, k). 



In view of the symmetries (2.3), this implies that ifi(x, t, k) and A<fi(x, t, k)A 1 satisfy 
the same ODE. Thus, the matrix J(k) defined by 

= e-^ k)x - i W t (Mi(x,t,k)- 1 AM 5 (x,t,uk)A~ 1 ) , keDi, (2.20) 

is independent of x and t. We have to show that J(k) = I. By (2.10), the matrices 
7 1 and 7 5 defined by ("J 1 )^ := 7^ and (7 )y := jfj are given by 
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Consequently, 
Mi(0,T,jfe) 




M 5 (0,T,uk) 




k G Di, 



(2.21) 



(2.22) 



where * denotes an unspecified entry. Evaluating (2.20) at (x,t) = (0, T) and using 



form 



(|2.22|) as well as the determinant condition detMi 

J(k) = 



1, we find that J(k) has the 



Similarly, evaluating (|2.20|) 
lim Mi(x,0,/c) = 



* * * 

* ■!= 
* != , 

00 and using that 




lim M$(x, 0,Lok) 



'1 * 
1 



k G Di 



,0 



1 
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we find that J(k) has the form 



J(k) 



1 N 

* 1 I . 

* * 0, 



It therefore only remains to show that J%2{k) = 0. Equation (2.21) yields 

(Mi(0,0,fc))32 = 0, (M B (0,0,wA;))i3 = 0. (2.23) 
Evaluating the relation 

- I 1 °\ 

Mi(x,t,A;)e £(fc):c+2(fe) * 1 = AM 5 {x, t, u}k)A~\ k € D u 

\0 J 32 (k) 1/ 

at (x,t) = (0,0) and using ( |2.23 ), we infer from the (32)-entry that 

J 3 2(fc)(Mi(0,0 J fc)) 3 3 = I keDt. 

Since J^ 2 {k) is an analytic function of k G D\ and (Mi(0, 0, k))^ — > 1 as k — > oo in 
Di, we conclude that J32 vanishes identically in D\. This proves (2.18) for k G D\. 
Equation (2.19) follows from the initial conditions 

Hj(xj,tj,k) = I, j = 1,2,3, 



and the symmetry properties (2.3) of L and Z in a similar, but simpler, way. □ 

The symmetry (2.18) implies in particular that 

S n (k) = AS n+4 (u}k)A~\ k G D n , n = 1,...,12, (2.24) 



and 



J n+1>n (x, t, k) = AJ n +B,n+4:(x, t, uk))A 1 , k G D n n D n+1 , n = 1 



where the index n is defined modulo 12. 



, . . . , 



12, 



(2.25) 



2.7. A matrix factorization problem. Define the 3 x 3-matrix valued spectral 
functions s(k) and S(k) by 

to(x,t,k) = fi 2 (x,t,k)e^ k)x+ ^ k)t s( y k), k G (ujS,oj 2 S,S), (2.26) 

m(x,t,k) = ii 2 (x,t,k)e^ k)x+i{k)t S(k), k G C, 



i.e. 



s(fc) = //3(0,0,AO, S , (fc)=/ii(0,0,fc). 



(2.27) 



We deduce from the properties of the /Zj's that s(/c) and S'(fc) have the following 
boundedness properties: 

s(k) is bounded and analytic for k G (ujS,uj 2 S,S), 

S(k) is bounded and analytic for k G {ujIZ U w7t, uj 2 1Z U o; 2 7£, 7£ U it), 

s (k) is bounded and analytic for k G (u>S,u) 2 S,S), 

S A (k) is bounded and analytic for k G (wT U uf, w 2 T U u; 2 7", Tut). 



Moreover, 



M n (x,t,k) = f i 2 (x,t,k)e^ x+ ^ t S n (k), k G D n . 



(2.28) 
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Proposition 2.6. The S n 's defined in (2.12) can be expressed in terms of the entries 
of s{k) and S(k) as follows: 



Si(k) 



S 3 (k) 



/Sli 
' Wi 

S 2 1 
W x 
S31 
Wi 



m 21 (s) 

S33 
mn(s) 





/ rn 22 (s) m 21 (S) 
m 12 (s) mnjS) 



Sl3\ 

«23 
S33/ 



S 2 (k) 



/ m 22 (S) m 2 i(s) 
Wo s 33 

m 12 (S) mn(s) 
Wi S33 




V 



S33 




W3 





S13 

•S23 
S33/ 



( m 22 (s) 5i2 

' 833 N Wl 

mi2(s) S 22 

S33 W4, 

u W4. 




(2.29a) 



(2.29b) 



where the functions {Wj(k)}f are defined by 



Wx(k) = S n m n (s) - S 21 m 21 (s) + S 31 m 31 (s) = (S T s A ) 



W 2 (k) = mu(s)m 22 (S) - m 2 x(s)m 12 (S), 
W 3 {k) = m 22 (s)mn(S) - m 12 (s)m 21 (S), 
Wi{k) = -S l2 m 12 (s) + S 22 m 22 {s) - S 32 m 32 {s) = (S T s A ) 22 . 
The functions {S n (k)}l 2 can be obtained from (2.29) and the symmetries in (2.24). 

Proof. It will be convenient to work with eigenfunctions which are defined for all val- 
ues of k G C. Thus, for each Xq > 0, let 7^° denote the contour (Xq, 0) — > (x, t) in the 
(x, i)-plane. We introduce fi 3 (x, t, k; Xq) as the solution of (2.15 ) with j = 3 and with 
the contour 73 repl aced by 73^°. Similarly, we define M n (x, t, k; Xo), n = 1, . . . , 12, 
as the solutions of (2.9) with 73 replaced by 73^°. We will first derive expressions for 
S n (k- X Q ) := M n (0,0,k;X ) in terms of S(k) and s(k; X ) := fi 3 (0,0, k; X ). Then 
(2.29) will follow by taking the limit Xq — > 00. 
Define R n (k; Xq) and T n (k; Xq), n = 1, . . . , 12, by 

R n (k; X ) = e- iT M„(0, T, k; X ), 

T n (k; Xo) = e-^ Xo M n (X , 0, k; X Q ), 

Then 

'M n {x, t, k; X ) = m(x, t, k)e tx+ ^R n (k; X ), 
< M n (x,t,k;X ) = ^ 2 (x,t,k)e^ +£t S n (k;Xo), 
M n (x, t, k; Xq) = n 3 (x, t, k; Xo)e tx+±t T n {k- X ), 



k G D n , 
k G D ri 

n = 1 



(2.30a) 
(2.30b) 

12, k G D n . 

(2.31) 



The relations (2.31) imply that 

s(k;X ) =S n (k;X )T- 1 (k;X ), S(k) = S n (k; X )R-\k; X ), k G D n . 

(2.32) 

These equations constitute a matrix factorization problem which, given s and S, can 
be solved for {R n , S n ,T n }. Indeed, the integral equations (2.9) together with the 
definitions of {R n , S n , T n } imply that 

(R n (k;X )) ij =Q if 75=71, 



(Sn(fc;*0))y 
(T n (k;X )) tJ 



if 7^=72 



(2.33) 



if ll 



73, 



where 7^ is defined by (2.10). It follows that (2.32) are 18 scalar equations for 18 
unknowns. By computing the explicit solution of this algebraic system, we find that 
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{S n (k; Xo)}f are given by the equation obtained from ( 2.29[ ) by replacing \S n (k), s (k)} 



with {S n (k; Xq), s(k; Xq)} . Taking Xq — > oo in this equation^] we arrive at (2.29). □ 



2.8. The global relation. The spectral functions S(k) and s(k) are not independent 



but satisfy an important relation. Indeed, it follows from (2.26) that 

/tii<r T S s = fj, 3 . 
Since fii(0, T, k) = I, evaluation at (0, T) yields the following global relation: 

5 _1 (fc)s(ife) = e-^ k)T c(T,k), k e (ojS,uj 2 S,S), (2.34) 

where the first, second, and thirds column vectors of the 3 x 3-matrix valued function 
c(T, k) := /X3(0, T, k) are analytic and of order 0(1/ k) as k — > oo for k in luS, uj 2 S, 
and 5, respectively. 

Remark 2.7. We also have 

M3 4 = A^e-^V 1 , ^ = tfe-e-tos*, 

and so 

^e-^-^C^-V) = M3 • 
Evaluating this equation at (0, T) we find the following cofactor version of the global 
relation: 

S T (k)s A (k) = e^ k)T c A (T,k), k G (ljS,u) 2 S,S), (2.35) 

where the first, second, and thirds column vectors of the 3 x 3-matrix valued function 
c A (T, k) := fi A (0, T, k) are analytic and of order 0(1/ k) as k — > oo for k in uaS, u 2 S, 
and 5, respectively. 



2.9. The residue conditions. Since 1x2 is an entire function, it follows from (2.28) 



that M can only have singularities at the points where the 5 n 's have singularities. 



In view of the symmetries of Lemma 2.5, it is enough to study the singularities for 



fc£5. We infer from the explicit formulas (2.29) that the possible singularities of M 
in S are as follows: 

• [M]2 could have poles in D\ U D2 at the zeros of S3s(k). 

• [M]i could have poles in D\ at the zeros of W\(k). 

• [M]i could have poles in D2 at the zeros of W2(k). 

• [M]i could have poles in D3 U D4 at the zeros of sss(k). 

• [M]2 could have poles in D3 at the zeros of W%(k). 

• [M]2 could have poles in D4 at the zeros of W^(k). 

We denote the above possible zeros by {kj}i and assume they satisfy the following 
assumption. 

Assumption 2.8. We assume that 

• s 33(&) has ni > simple zeros in D\ U D2 denoted by {/cj}™ 1 , 

• W\(k) has ri2 — n\ > simple zeros in D\ denoted by 

• Wi(k) has 713 — ri2 > simple zeros in D2 denoted by {A;j}^ +1 , 

• sss(k) has — 713 > simple zeros in D3 U D4 denoted by {kj} T n3+1 , 

• Ws(k) has n5 — m > simple zeros in D3 denoted by {kj}™^ +1 , 



1 T14 



W^(k) has — 715 > simple zeros in D4 denoted by {kj}^ 



+i' 



9 

Note that all quantities in this equation have well-defined limits as Xq — > 00. Had we instead 



tried to let Xq — > 00 already in ( 2.32 1 , the limit of the first equation in (2.321 would only have been 
well-defined for k € (uS, uj 2 S, S). This is the reason for introducing Xq- 
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and that none of these zeros coincide. Moreover, we assume that none of these 
functions have zeros on the boundaries of the -D n 's. 

In the next proposition we determine the residue conditions at these zeros. 



Proposition 2.9. Let {M n }\ 2 be the eigenf unctions defin ed b y (2.9) and assume 
that the set {kj}^ of singularities in S are as in assumption 2.8. Then the following 
residue conditions hold: 



Res [M]r 



Res [M]i 



m 2l {s{k j )) g 32fe 



e^W'lMikj)^, 1 < j < n u kj GDiU D 2 , (2.36a) 



W 1 {k j 



1 pll(fcj>23(fcj) ~ ' S 2l(M Sl ^M e fel(fcj)[ M ( fcj .)] 2 



m 3 i(s(kj)) 



^e^^)[M(fc i )] 3 }, m<j< n 2 , k 3 G A 



S33(%) 

^' M ]'= rl?' )) rf» e ' ,2 ''' iji ' M ft»^ 



(2.36b) 

n 2 < j < n 3 , fcj G D 2 , (2.36c) 



Res t M ]i = 7f ( v^l\ ee3lfe) [ M (^)]3, n 3 <j< n 4 , % GD 3 U D 4 , (2.36d) 
kj s 13 (kj)sz 3 (kj) 



Res [M] 2 



w^kj^Mkj)) 
+ 



n 4 < j < n 5 , kj E D 3 , (2.36e) 



1 f Si 2 {kj)s 23 {kj) - S 22 {kj)si 3 {kj) cdi2{ki) k 

• ^\ m 32 {s{kj)) ' 3 1 

s 33 (kj) 



e»Wi)[M(kj)] 3 }, n 5 <j<N, kj G D A 



(2.36f) 



where f := df jdk, 9ij is defined by 

Oij(x, t, k) = (k{k) - lj(k))x + (z t (k) - Zj(k))t, i,j = 1, 2, 3, (2.37) 
and we have suppressed the (x,t) -dependence for simplicity. 



Proof. We will prove ( 2.36a )-( 2.36c); the other conditions follow by similar argu- 
ments. Equation (2.28) implies the relations 

Mi = ^ 2 e ix+it Si and M 2 = ^ 2 e ix+ ^S 2 . (2.38) 

In view of the expressions for S\ and S 2 given in (2.29a), the three columns of ( 2.38^ ) 
read 



[MJ 2 e»" Ml + =H%2] 2 , 

= si 3 e 9l3 [/i 2 ]i + s 23 e e ^ 2 } 2 + s 33 [^ 2 } 3 , 
while the three columns of ( 2.38b ) read 



2j3 ; 



[M, 



2 1 



m 22 (S) 



LA*2jl H 777 — e 21 [/"2j2 



W 2 Lr ^ J± M^ 2 
[M 2 ] 2 = m21 Me^Ni + mil(s) 



[^2, 



S33 S33 

[M 2 ] 3 = Sise 913 ^]! + s 23 e 023 [/i 2 ] 2 + S33H3. 



(2.39a) 

(2.39b) 
(2.39c) 

(2.40a) 

(2.40b) 
(2.40c) 
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In order to prove (2.36a), we first suppose that kj G D\ is a simple zero of s 33 (k). 
Solving (2.39c) for [/^ji and substituting the result into ( 2.39b| ), we find 



[M- 



12 



S13S33 



S13 



1 m 31 {s) 

A*2 3 [fJ-2 2- 

S13 



Taking the residue of this equation at kj , we find the condition (2.36a) in the case 
when kj G D\. Similarly, solving (2.40c), substituting the result into (2.40b), and 
taking the residue at kj, we find that (2.36a) also holds if kj is a simple zero of S33 
in D 2 . 

In order to prove (2.36b), we suppose that kj G D\ is a simple zero of Wi(A;). 
Solving (2.39b) and (2.39c) for [^h and [^2)2 and substituting the result into ( 2.39a| , 
we find 



[M 



0831 



11 



'S'll523 — £21313 Q 21 r 



m 3 i(s) IVim 3 i(s) \^l s 33 S33"l 3 i(Sj 



£31 



e S31 [Mi] 3 . 



(2.41) 



Taking the residue of this equation at kj, we find (2.36b). 

In order to prove (2.36c), we suppose that kj G D 2 is a simple zero of W2(&)- 
Solving (2.40b) for [/J.2J 1 and substituting the result into (2.40a), we find 

P e 2 i 



[M, 



2 1 



7tN2 + TT r T~\ M 2 2 

m 2 i{s) W 2 m 21 (s) 



(2.42) 



Taking the residue of this equation at kj, we find (2.36c). 



□ 



3. The Riemann-Hilbert problem 



The sectionally analytic function M(x, t, k) defined in section[2]satisfies a Riemann- 
Hilbert problem which can be formulated in terms of the initial and boundary values of 
q(x, t) and r(x, t). By solving this RH problem, the solution of ( |1.1[ ) can be recovered 
for all values of x, t. 



Theorem 3.1. Suppose that {q(x,t),r(x,t)} is a solution of the system (1.1) in the 



half-line domain {0 < x < 00, < t < T} with sufficient smoothness and decay 
as x — )• 00. Then u(x,t) can be reconstructed from the initial and boundary values 



{q (x),r (x),g (t),h (t),g 1 (t),hi(t)} defined in (1.4) as follows. 



Use the initial and boundary data to define the jump matrices J m)H {x,t, k), n,m = 
1, . . . , 12, by equation (2.14) as well as the spectral functions s(k) and S(k) by equation 



(2.21). Assume that the possible zeros {kj}i of the functions S3 3 (/c) and {Wj(k)}f 



in S are as in assumption \K_ 

Then the solution {q(x,t),r(x,t)} is given by 



q{x,t) 



iV3 lim (kM(x,t,k)) 2 3, 

fc— >OD 



r(x, t) 



-VjJ' 



'V3 lim (kM(x,t,k)) 13 , (3.1) 



where M(x, t, k) satisfies the following 3x3 matrix RH problem: 

• M is sectionally meromorphic on the Riemann k-sphere with jumps across 
the contours D n n D m , n, m = 1, . . . , 12, see Figure^ 



Across the contours D n nD T/ 



n, m 



12, M satisfies the jump condition 



(2A§. 

M(x,t,k) 



I + O 



00. 
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The first column of M has simple poles at k = kj for m < j < n±. The second 
column of M has simple poles at k = kj for 1 < j < n\ and < j < N. The 



associated residues satisfy the relations in (2.36). 
For each zero kj in S, there are two additional points, 



ujk 



0J 2 k 



3 1 



at which M also has simple poles. The associated residues satisfy the residue 



conditions obtained from (2.36) and the symmetries (2.18). 



Proof. It only remains to prove (3.1) and this equation follows from the large k 
asymptotics of the eigenfunctions, see subsection 5.2 below for details. □ 



4. LlNEARIZABLE BOUNDARY CONDITIONS 



Theorem |3. 1| expresses the solution {q(x, t), r(x, t)} of the system (1.1) in terms of 
the solution of a RH problem whose jump matrix is given in terms of the spectral 
functions s(k) and S(k). The function s(k) is defined in terms of the initial data 
{qo( x ), r o( x )} via a system of linear Volterra integral equations, whereas the function 
S(k) is defined in terms of the boundary data {go(t), ho(t), gi(t), hi(t)} also via a 
system of linear Volterra integral equations. However, for a well-posed problem, only 
part of the boundary data can be independently prescribed; for example, for the 
Dirichlet problem go an d ho are prescribed, whereas for the Neumann problem g\ 
and hi are prescribed. 

In the next section, we will show that in general S(k) can be determined from 
the given boundary conditions and the initial data via a system of nonlinear integral 
equations. In this section, we will consider the special case of linearizable boundary 
conditions. These are boundary conditions for which it is possible to express S(k) in 
terms of s(k) and the given boundary conditions via algebraic manipulations. Thus, 
the above nonlinear step can be avoided and the problem can be fully solved by linear 
operations alone. 

Our present goal is to find an approach for analyzing linearizable boundary condi- 
tions for equations with 3x3 Lax pairs. To this end, we will consider IBV problems 



for the system (1.1) with arbitrary initial data and vanishing Dirichlet boundary 
conditions: 

9(0, t) = r(0, t) = 0, < t < oo. (4.1) 

We will show that these particular boundary conditions are linearizable in the sense 
that the jump matrix as well as the residue conditions of the RH problem of Theorem 



3.1 can be expressed entirely in terms of s(k). The elimination of S(k) is achieved 
by analyzing the global relation and the equations obtained from the global relation 
under those transformations in the complex /c-plane which leave the associated dis- 



persion relation invariant. For the system (1.1), there is one such transformation, 
namely k i— > —k (which leaves k 2 invariant). 

For a function f(k), we define f(k) by f(k) = f(—k). 



Theorem 4.1. Let {q(x,t),r(x,t)} satisfy the system (1.1) in the half-line domain 
{0 < x < oo, < t < oo} together with the initial conditions 

q(x, 0) = qo(x), r(x, 0) = ro(x), < x < oo, 

and the Dirichlet boundary conditions 

q(0, t) = r(0, t) = 0, < t < oo. (4.2) 
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Assume that the initial and boundary conditions are compatible at (x,t) = (0,0). 



Then the jump matrices J m ,n{k) for the RH problem in Theorem 
explicitly in terms of the spectral function s(k) as follows: 



3.1 



are given 



J 



1,2 



J: 



3,4 



1 o o\ 

10 

e 9 3is 3 i n | 
S33W1 U ' ' 



1 

n e g 32g 32 -. 

S33W4 



J: 



2,3 



V 



1 

e 6 2l(s T s A )i2 
Wi 




(S T s A ) 21 (s T s A ) 12 
W1W4 





4,5 



/ 1 



V 



S33 




1 







ra 22 (s) 
S11S33 



/ 



(4.3) 



where 9ij(x,t,k) is given by (2.31) and the functions {Wi(A;), W 4 (&)} are defined by 
Wi (fc) = snmii (s) - s 2 im 2 i (s) + s 3 ira 3 i (s) , (4.4a) 
W 4 (/c) = -Sl 2 mi 2 (s) + S22«1-22(s) ~ ■532m 32 (s). (4.4b) 

XTie expressions for the remaining jump matrices follow from the expressions in U- 



and the symmetry (2.25). Moreover, the residue conditions (2.36) are given in terms 
of s(k) by 



Res [M] 2 



Res [M]t 

kj 



Res [M]i 



Res [M]i 

kj 



Res [M] 2 



Res [M] 2 



m 2 i(s(k j )) g 32(fc . 



Sl3(fci)»33(%) 



e e32 ^^[M(A; j )]3, 1 < j < m, kj GDjU D 2 , (4.5a) 



1 



+ 



Sll(fcj)g23(fcj) ~ S2l(fcj)gl3(fcj) fl 21 (fc. 



m 3l(s(^')) 



[M(^-)]2 



§ ^ e S 31 (k j)[M{kj)h 



n\ < j < n 2 , kj G Z>i, 
• ; /(%)]2, 



S33(fcj) 

sii(fcj)a 3 3(fcj) ~ a 3 i(A;j)ai3(fcj) ^ 21 (fc . } [ 

VVi(fcj)»n2i(«(A;j)) 

n 2 < j < n 3 , kj G D 2 , 

m 22 (g(fcj)) e 3l(fc .) 



(4.5b) 



(4.5c) 



^3i^J[ M(A ..)] 3) n 3 < j < n 4 , fcj G D3 U D 4 , (4.5d) 



Sl3(^i)s33(^i) 

Sl2(fej)s 33 (fcj) - a32(fcj)ai3(fcj) p 12 (fc 

W4(A;j)m 22 (s(/c i )) 



c- "vj)[M(kj)]x, 

m < j < TI5, kj G D 3 , 
[M(kj)] x 



+ s 32(fci) c e 32(fc ,) 

S33(fcj) 



gi2(%)s23(fcj) ~ a22(fcj)ai3(fcj) e e 12 (fc . ) 
m 32 (s(*y)) 

[M(^)] 3 ), 



"5 < j < fc,- G -D4. 



(4.5e) 



(4.5f) 



Proof. The vanishing boundary conditions (4.2) imply that the matrix Z defined in 
(2.1 ) satisfies Z(x, t, k) = Z(x, t, —k). It follows that the spectral function S(k) obeys 
the symmetry 



S(k) = S{-k), 



k G C. 



(4.6) 



We will use this symmetry property together with the global relation (2.34) to elim- 
inate all the entries of S(k) from the jump matrices and the residue conditions. 

Since we are now assuming that T = 00, the columns of S(k) have smaller domains 
of definition than above (but the domains of boundedness remain the same). More 
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precisely, 

S(k) is defined and bounded for k G (ujIZ U uTZ, uj 2 1Z U uj 2 TZ, 1ZU71). 
The global relation when T = oo is 

, A; G 
C33/ 



(5 




L> 7 UD 8 
D 5 UD 6 



D 9 U L>i 


D n U D12 



L> 3 U D4 
D l VJD 2 



(4.7) 



(4.8) 



where 



(S A ) T (k) 



mn(5) 
-mi 2 (5) 
77113(5*) 



-77721 (S) 

m 22 (5) 
-m 23 (5) 



m 31 (S) 
-m 32 (5) 



and the notation indicates that the (11) entry of (4.8) is not valid for any k G C, the 



(12) entry of (4.8) is valid for k G D g U Di etc. 



We first derive the expression for J 2j 3. From the explicit expressions in (2.29), we 
see that the matrices S 2 (k) and S 3 (k) only depend on the entries of S(k) via the 
quotients 

m 12 (S) mn (5) 



m 22 (S) 



and 



77721 (S) ' 



yields two 



respectively. Because of the symmetry (4.6), the global relation (4. 
algebraic relations for the entries of S(k) in each of the Dj 7 s. In D 2 , these relations 
are 



-si3mi 2 (5) + S23m 2 2(5') 
-sum 12 (5) + S2im 2 2(5') 



S33m 32 (5) 
S3im 32 (5) 



0. 
0. 



(4.9a) 
(4.9b) 



Equation (4.9a) is the (23) entry of (4.8 ), whereas (4.9b) is obtained by letting k — > —k 
in the (21) entry of (4.8) and using (4.6) in the resulting equation. The two relations 
in (4.9) can be solved for the two quotients 



mi 2 (S) 
m 32 (5) 



and 



m 22 (S) 
m 32 (5) 



with the result that 



mi 2 (5) 
m 32 (5) 
m 22 (S) 



S23S31 - S21S33 



S13S21 
S13S31 



S11S23 
S11S33 



k G D 2 , 
k G D 2 . 



(4.10a) 
(4.10b) 



m 32 (S) S13S21 - sns 23 
Using these relations to eliminate the S(k) dependence in the expression for S 2 in 



(2.29a), we find 



S 2 (k) 



S11S33-S13S31 

S33 (SlH7lll(s)-S2lW2l(«)+S3lW3l(s)) 

S21S33-S23^31 

S33(siH7lll(s)-S2lW2l(«)+S3lW3l(s)) 





(4.11) 



Similarly, the (13) entry of (4.8) together with the equation obtained by letting 
—> —k in the (12) entry of (4.8) yield two relations for the entries of S(k) in .D3. 



k 

Solving these two relations for the two quotients 

mn (5) , m 2 i(5) 
and 



77731(5) 



77731(5)' 
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we find 



m n (S) _ S23S32 - S22S33 
m 3 i(S) S13S22 - S12S23' 
m,2i(S) _ S13S32 - S12S33 
m 3 i(S) S13S22 - si2S23 : 



k G D 3 



(4.12a) 
(4.12b) 



Using these relations to eliminate the S(k) dependence in the expression for S3 in 



(2.29b), we find 



S 3 (k) 



/ m 2 2(s) 
m 12 {s) 



V 



S33 




S33(3l2mi2(s)-.S22m22(s)+S32«l32(s)) 

«23^32— ^22^33 

■S33(*12mi2(s)-S22m22(s)+S 3 2m 3 2(s)) 





(4.13) 



he expression for J23 in (4.3) now follows by substituting the matrices (4.11) and 



(4.13) into the definition J2; 



Cx+Zt 



(iS 2 ~ 1 S3), recalling that both matrices have 
unit determinant. 

We next derive the expression for J\ %. The approach used for finding ^3 is not 
successful in this case (see remark 4.2 below), and we have to rely on a somewhat 
different argument. 



The explicit expression for the matrix 82(h) in (2.29a) only involves entries from 
the following column vectors: 

Ms, [s A h, [s%. 

Since all these column vectors are analytic in D\, 82(h) admits an analytic continu- 
ation to D\. It follows that the defining equation 

J 12 = e^ +it (S^S2), (4.14) 

which a priori is valid only for h G Di(lD2, can be extended to the line h £ D\2^D\. 



Substituting the explicit expressions in (2.29) for S\(h) and 82(h) into (4.14), we find 



after some algebra that depends on the entries of S(k) only via the four quotients 
m 13 (S) m 2 3(S) mi 2 (S) m 2 2(S) 



ra 33 (S)' rri33(S)' m 32 (5)' m 32 (S) 



(4.15) 



The (32) entry of (4.8) together with the equation obtained by letting h — > —h in 



the (31) entry of (4.8) yield two relations for the entries of S(h) in D\2- Solving these 



two relations for the first pair of quotients in (4.15), we find 

mi3(S) S22S31 - S21S32 



m-33(S) 
m 2 3(S) 



S12S21 
S12S31 



S11S22 
S11S32 



k € D 



12, 



k G D 



12- 



(4.16) 
(4.17) 



m 3 3(S) S12S21 - S11S22' 

We use these quotients to eliminate the first two quotients in (4.15) from the expres- 
sion for Ji t 2 for k G D12 H D\. 

On the other hand, since all quantities in the equations (4.10) are analytic in D\ 



the relations in (4.10) are valid also for k G D\% n D± and can be used to eliminate 



the second pair of quotients in (4.15) from the expression for for h G Z?i 2 D D\. 



This yields after simplification the following expression for for h G Dyi n D\\ 

1 ON 

1 

1, 



h G D\2 Pi D\. 



Lg31 



833(ail"lll(s)-S2lTO2i(s)+83im3i(s)) 
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Since all quantities on the right-hand side are analytic in D\ , this expression for J\ 2 



is valid also for k G D\ n D 2 . This establishes the expression for J± j2 in (4.3) 



The proof of the expression for J 3i 4 is similar: Since the expression (2.29b) for 
S 3 {k) is analytic in D4, the defining relation 



J: 



3,4 



Xx+Zt 



(S 3 1 S 4 ), 



which a priori is valid only for k G D 3 n D4, can be extended to the line k G D4 n D^. 
After some algebra we find that depends on the entries of S(k) only via the four 
quotients 

m 13 (S) m 23 (S) m u (S) m 21 (S) 



m 33 (S)' m 33 (S')' m 3 i(S)'' m 3 i(S)' 



Solving the (31) entry of (4.8) and the equation obtained by letting k 
(32) entry of (|4.8[) for the first two of these quotients, we find 



mi 3 (S) _ S22S31 - s 2 is 32 
m 33 (S) S12S21 - s n s 22 ' 
m 23 (S) _ S12S31 ~ sng32 
"133(5) 



S12S21 - S11S22 



k G D, 



k G D 5 . 



(4.18) 
k in the 

(4.19a) 
(4.19b) 



We use these relations to eliminate the first pair of quotients in (4.18) from the 



expression for J 3) 4 for k G -D4 n -D5. We use the relations in (4.12), which admit an 



analytic continuation to D4, to eliminate the second pair of quotients in (4.18). This 



yields the expression for J 3j 4 given in (4.3), but for k G -D4 D -D5. However, since the 



expression for J 3j 4 in (4.3) is analytic in D4, the equation is also valid for k G -D 3 n.D4. 



Finally, the expression for follows simply by substituting the expressions for 



S4 and 55 in (2.29) into the definition 



Xx+Zt 



(5 4 1 S5); 



in this case the S(k) dependence disappears automatically. 



We now show that the residue conditions in (2.36) can be written as in (4.5). 



In order to show (4.5e), we note that 
m 2 i(S(kj)) m 2 i(S) 



W 3 {k,j 



Res 



W 3 



k mii(g) 



ni2 



^m 22 (s) - mi 2 {s) 



where kj G D 3 is the zero of W 3 . Utilizing (4.12) to eliminate the S(k) dependence 



from this expression, we find after simplification 



m 2 i(S(kj)) 



Res 



S12S33 ~ ^13^32 
S 33 W4 



where W4 is given by (4.4b). Substitution of this into (2.36e) yields (4.5e). 



In order to show (4.5f), we use equations (4.12) and (4.19). As we have seen, both 



of these equations are valid for k G D4 n D5. Algebraic manipulation using these 
relations shows that 



'S'l2'S23 — Si 3 S 22 S12S23 — «13^22 



m 32 {s)W 4: 



m 3 2{s)W4 



k G D 4 HD 5 , 



and 



'32 



S32 



S33W4 S33W4 ' 



k G -D4 n D 5 . 
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Since all quantities in these equations are analytic in D4, these relations hold for all 
k G D4. In particular, they hold at kj G D4. In view of (2.36f), this yields (4.5f). 
The proofs of (4.5b) and (4.5c) are similar. □ 



Remark 4.2. It is tempting to attempt to derive the expressions for the matrices 
Ji,2 and J 3 ^ in Theorem 4.1 in the same way that the expression for J 2}3 was derived. 
However, this does not appear to be possible. Indeed, consider t he de rivation of 
the expression for J\p = e Cx+zt (S^ 1 S 2 ) . In view of the expression (4.11) for S 2 (k), 



which is independent of S(k), we only need to find an analogous expression for S\(k). 



Prom (2.29a) we see that Si(fc) only depends on S(k) via the quotients Su/S 3 i and 
S21/S31. We can use the cofactor version of the global relation, which for T = 00 is 
given by 



S T s A 





D 7 

I>io 
D 8 U A), 

(4.20) 



to eliminate these quotients from Si(k). The (11) entry of (4.20) and the equation 
obtained by letting k — > —k in the (13) entry of (4.20) provide two algebraic relations 
which can be solved for Sn/S 3 i and S21/S31 in D\. The result is 



S u -m 31 (s)m 23 (s) + m 21 (s)m 33 (s) + ^^-m 23 (s 



S2I 



mii(s)m 23 (s) - m 2 i(s)mi 3 (s) 
■m 31 (s)m 13 (s) + m u (s)m 33 (s) + m || c) mi 3 (s) 



k G D 1 , 



k£D 



mn(s)m 23 (s) - m2i(s)mi 3 (s) 

We can use these relations to eliminate S11/S31 and S21/S31 from Si(k). However, the 
result still involves the unknown quantity ^5^, hence the solution is not effective. 



5. NON-LINEARIZABLE BOUNDARY CONDITIONS 

A major difficulty of initial-boundary value problems is that some of the boundary 
values are unknown for a well-posed problem. On the other hand, all boundary 
values are needed for the definition of S(k), and hence for the formulation of the 
RH problem. In the previous section we analyzed the special case of linearizable 
boundary conditions — these are boundary conditions for which the matrix S(k) can 
be eliminated algebraically from the formulation of the RH problem. In general, 
for non-linearizable boundary conditions, we need to determine S(k) from the initial 
data and the given boundary values. In this section we concentrate on the effective 
characterization of S(k) . Following [T7] , we call a characterization of S(k) effective if it 
fulfills the following requirements: (i) In the linear limit, it yields an effective solution 
of the linearized boundary value problem, (ii) For 'small' boundary conditions, it 
yields an effective perturbative scheme, i.e. it yields an expression in which each 
term can be computed uniquely in a well-defined recursive scheme. 

Our main result expresses the spectral function S(k) in terms of the prescribed 
boundary data and the initial data via the solution of a system of nonlinear integral 
equations. As in [T7], our approach uses three ingredients: (a) The large k asymptotics 
of the eigenfunctions. (b) The global relation, (c) A perturbative scheme to establish 
effectiveness. 
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5.1. The global relation. Evaluating (2.26) at (x,t) = (0,t), we find 

H 2 (0,t,k)e zt s(k) = c(t,k), k G (uS, u 2 S, S), 
where c(t,k) = fi 3 (0,t,k). Denning functions {$>j(t, k)}f and {cj(t,k)}\ by 



fi 2 (0,t,k) 



'<S> 3 (t,LJ 2 k) <S> 2 (t,uk) $l(t,k) S 

$i(i,u; 2 A;) $ 3 (t,u)k) $ 2 (t,k) 
$ 2 (t,uj 2 k) $i(i,wfc) $ 3 (f,A:), 



[c(t, fc)k- 

S33(^) 



c 2 (t,k) I . 
K C3(t,k) 



we can write the (13) and (23) entries of this equation as 

<Z> 3 (t,oj 2 k)e {zi - Z3)t — + & 2 {t, cok)e {z2 - Z3)t — + $ x (i, fc) = ci(t, fc), 



«33 



S33 



w 2 *)^* 1 -* 3 ^ + $ 3 (t, wfcje^-* 3 ^^ + $ 2 (*, fc) = c 2 (i, fe), 
•S33 S33 



fc G S, 
(5.1a) 

fc G 5. 
(5.1b) 

The functions {cj(t, fc)} 3 are analytic and bounded in S away from the possible zeros 
of S3 3 (k) and of order 0(1/ k) as k — >■ 00 in 5. The functi ons fc)}f are entire 

functions of which are bounded for k G 7~U7". Equation (5.1a) shows that whereas 
each of the functions $ 3 (i, u} 2 k)e^ Zl ~ Z3 ^, ^ 2 (t,ojk)e^ Z2 ~ z ^ t , an d <&i (£, fc) is bounded 
in TUT, the combination appearing on the left-hand side of (5.1a) is bounded also 
in 5 (away from the zeros of s 33 ). A similar remark applies to (5.1b). 



5.2. Asymptotics. An analysis of (2.5) shows that the eigenfunctions {/ij}f have 
the following asymptotics as k — > 00: 



fij(x,t,k) = I + 



(x,t) 



{Xj,tj) 



A J 



+ 



(Xj ,tj ) 



1 



urr ujq 

? 2 ,.,,2 



r 2 oj 2 q 2 




(5.2) 








1 







,U 2 Q X 




— n(— 

k 2 + U 3 



A; 



00, 



3 = 1,2,3, 



where each column is valid within its region of boundedness and the closed one-forms 
A and {%}? are defined by 



A = qrdx + 



i (q x r - qr x ) 2 , 3 3 , 

— 71 3< fl +r) 



(*,*) x 2 
A 

(Xj ,tj ) 

q 2 r 2 + -g^ - \(rq xx + qr xx 



6 



i (Qxr - gr x ) 
2^/3 



dt, 



(q 3 + r 3 ) 



1,2,3. 



dx 



Hence 

[s(k)h 







, .2 r(0,0) a 



i\/3 



W(0,(Vf 
-9(0,0) 




(5.3) 
00, k G S, 
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and 



*l(t,fc) 



k k 2 ^ u 3 



(fc 3 ) ! 



fc -)• oo, k G TUT, 

fc — >■ oo, fe e TuT, 
fc — >■ oo, fe e TuT, 



(5.4a) 
(5.4b) 
(5.4c) 



where 



go(t) 

»V3' 



*f°(t) 



3 + ^3 



A + 



(0,0) 



*?>(t) 



w/*,(t) 2 kAo /' (0 ' t) a , Si(t) 



a/5 



A + 



(0,0) 



(o,t) 



A. 



(0,0) 



In particular, we find the following expressions for the boundary values: 



5o = « 



2 ' 



-ltd 



gx = 3$[, 2) - tdhl + iV3g <S> { 3 1] , h ± = 3w$f } 
We will also need the asymptotics of ci and C2- 



(5.5a) 
(5.5b) 



Lemma 5.1. The global relation (5.1) implies that the functions c\ and c 2 satisfy 

— > oo, G 5, (5.6a) 



ci(tjt) ^ + £ ^) +0( : 



jfe 



— > oo, k G 5. 



(5.6b) 



Proof It follows from ( 2.5 ) that the $j 's admit an expansion of the form (see chapter 
6 of 0) 



* 2 (t, fc) = (ao(t) + -4 1 + • ■ • ) + (A) (*) + 



+ 



,(*l-«3)* 



(5.7) 



+ 7o(*) + 



71 0) 



k 



+ ■■■ e 



D (z 2 -z 3 )t 



k —> oo, k G 



where the coefficients aj(t), (3j(t), Jj(t), j = 0, 1, . . . , are column vectors which are 
independent of k. We determine the coefficients by substituting (5.7) into the system 



ktd 9o + — - h Z Q j $ 2 + U^O - 50 - -y| ) *3, 



$2t + (z3 - z 2 )^2 = ktdh - go 



= ( kd 2 g + - h 2 Q ) $i + ( kh - g 2 - ^ ) $ 2 , 



ih\ 



V3 



and using the initial conditions 



- ! % + ^|-/io I*- 



a (0)+/?o(0)+7o(0) = (0,0, 1) T , ai(0)+A(0) + 7i(0) = (0,0, 0) T . 
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This yields 



'*i(t,fc) 




(zi(k)-z 3 (k))t 



(z 2 {k)-z 3 (k))t 



k — > oo, 6 



Substituting this expansion and the expansion (5.3) of [s]s into (5.1a), the result- 
ing left-hand side involves the exponentials e*- 21-23 ^ and e^ 22-23 - 1 '. Since the global 
relation requires that the right-hand side is of order 0(l/k) as k — > oo in S, the 



coefficients of these exponentials must vanish. The remaining terms yield (5.6a). The 
expansion (5.6b) follows in a similar way from (5.1b). □ 



5.3. The Dirichlet and Neumann problems. The following theorem expresses 
the spectral functions A{k) and B{k) in terms of the prescribed boundary data and 
the initial data via the solution of a system of nonlinear integral equations. 

Theorem 5.2. Let T < oo. Let qo(x), x > 0, be a function of Schwartz class. For 
the Dirichlet problem it is assumed that the functions go(t) and ha(t), < t < T, 
have sufficient smoothness and are compatible with qo(x) at x = t = 0. Similarly, for 
the Neumann problem it is assumed that the functions g\{t) and hi(t), < t < T, 
have sufficient smoothness and are compatible with qo(x) at x = t = 0. Suppose that 



S33(^) is free from zeros (see remark 5.3 below for the case when s^{k) has zeros). 



Then the spectral function S{k) is given by 

(A{uj 2 k) B{ujk) C(k)\ 

S(k) = C(u 2 k) A(u;k) B(k) , (5.8) 
\B(oj 2 k) C(u>k) A(k) J 

where 

A(k) = * 3 (T, w 2 £;)$3(T, uk) - $i(T, u} 2 k)<S> 2 (T, uk), (5.9a) 

B(k) = - l$z{T,u 2 k)<S> 2 {T,k) - $i(T,u; 2 £;)$i(T,£;)] e i^)-z 2 (k))T ^ ^ ^ 

C{k) = [^(r.wfcJ^CT.fcJ-^aC^w^^iCT.fcJle^W-* 1 ^, (5.9c) 
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and the complex-valued functions {$j(t, k)}f satisfy the following system of integral 
equations: 



*l(t,fc) 



,( Zl (k)-z 3 (k))(t-t>) 



kujg + ^= - hi j $ 2 
+ ^ 2 /i - gg - ^0 



,(^(fc)-^3(fc))(t-t') 



few/to - g 



ih] 



$ 3 (t',k)dt', (5.10a) 



[t',k)dt\ (5.10b) 



* 3 (t ) fc) = l + 



(t',k)dt', 

(5.10c) 

For a function f(k), let f+(k) and f-(k) denote the following even and odd com- 
binations of f(k) : 

f + (k) = f(k) + f(-k), = f(k) - f(-k), k G C. 

(a) For i/ie Dirichlet problem, the unknown Neumann boundary values gi(t) and 
h\(t) are given by 

3V3g (t) 



9i(t) 



hi(t) 







"-/ 


k 


™ Jc 


s . 


2,00 r 

2vri /as 


k§ 


3u j 


k 

3S 


m J 



2g (t) 
iV3 



dk - uho(t) + 



2tt 



$3-(t,k)dk 



2ojh (t) 



$i(t, -w 2 /t)e^- 23 ^^ + * 3 (i, -wfcW 

S33 «33 



dk-u 2 g 2 (t) 



S 

3V3uj 2 h (t) 



2tt 

(z 2 -z 3 )t^3 



dk, (5.11a) 
<S> 3 ^(t,k)dk 



us 



$ 3 (t, -w-W* 1- * 3 ' 1 ^ + $ 2 (*, -wfe)e 

•S33 S33 



dfc. (5.11b) 



(b) For i/ie Neumann problem, the unknown Dirichlet boundary values go(t) and 
ho(t) are given by 



flo(t) 



2vr 



r)5 



$ 2 +(t,k)dk 



(5.12a) 



/io(t) 



+ 



^ ./as 
2vr 

V5 



$ x (i, -w 2 A;)e (2l - 23)i — + $ 3 (t, -wfc)e 
S33 



( a2 _ Z3 ) t S23 
•§33 



95 



$ 2 +(i, fc)dfc 



(5.12b) 



7T 



$ 3 (t, -w 2 fc)e (zi - 23) *^ + $ 2 (-wA:)e ( 

S33 S33 



= ( 22 - Z3 )if23 



dife. 



Proof. The representations (5.9) follow from the relation S(k) = e ZT fi^iOjT, k) T . 
(a) In order to derive (5.11) we note that equation (5.5b) expresses g\ in terms of 
and <3?2 2 ' '• Furthermore, equations (5.4) and Cauchy's theorem imply 



~<t> < i\t)= I [«3(t,fc)-i]dfc 



9T 



[$ 3 (t,fc)-l]dfc 



(5.13) 
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dT 



dk 



dT 



k$ 2 (t,k) - ^\t) 



dk. (5.14) 



Thus, 
2i?r 
~3~ 



and 



® < i\t) = -(l + [ )mt,k)-i]dk= ( [ + f )[$ 3 {t,k)-i\dk 

\JdT JafJ \Jas JasJ 

= [ [& 3 (t,k)-l]dk- I [& 3 (t,-k)-l]dk= [ <$>3-(t,k)dk (5.15) 
Jas Jas Jas 



2m 



*?>(*) 



+ 

as Jas 



[k$ 2 (t,k) - (t)]dk 



as Jas 



dk + I(t) 



as 



[k$2-(t,k)-2Q>£\t)]dk + I(t), 



(5.16) 



where I(t) is defined by 



I(t) = 2 [ [k$ 2 (t, k) - ^ (t)] dk. 
Jas 



The last step involves using the global relation (5.1b) to compute I(t): 
I(t) = 2 [ [kc 2 (t,k) - <5> { 2 \t)}dk 



as 



(5.17) 



dk 



k $i(tV*:)e (zi - 23) *— + $ 3 (t,^)e (z2 - 23) * — 

dS L s 33 S33. 

Hence, using Cauchy's theorem and the asymptotics ( 5.6b| ) of c 2 to compute the first 
term on the right-hand side of (5.17) and using the transformation k — > — k in the 
second term on the right-hand side of (5.17), we find 

4-7TZ 



I(t) 



US 



•S33 s 33 



dk. 



(5.18) 



Using (5.15), (5.16), and (5.18) in (5.5b) we find (5.11a). The proof of (5.11b) uses 
(5.1a) and is similar. 



(b) In order to derive (5.12) we note that equation (5.5a) expresses go in terms of 

(5.19) 



$2 • Furthermore, equations (5.4) and Cauchy's theorem imply 

$ 2 (t,k)dk 



dT 



dT 



$ 2 (t,k)dk. 



Thus, 



2m (i) 



+ 

aT Jar 



$ 2 (t,k)dk = 
$ 2 (t,k)dk + K(t) 



+ 

as Jas 



$ 2 (t,k)dk 



us 



as Jas, 
® 2+ (t,k)dk + K(t), 



(5.20) 
(5.21) 
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Re k 



Figure 4. The contour T defined in such a way that —T passes above the zeros of 
s 33(k) in S. The zeros of 833(h) are indicated by x's in the figure. 



where K(t) is defined by 



K(t) 



2 / <S> 2 (t,k)dk. 
'OS 



The last step involves using the global relation (5.1b) to compute K(t): 



K(t) = 2 c 2 (t,k)dk 
JdS 



(5.22) 



as 



dk. 



•S33 S33. 
Hence, using Cauchy's theorem and the asymptotics ( 5.6b| ) of c 2 to compute the first 
term on the right-hand side of (5.22) and using the transformation k — > — k in the 



second term on the right-hand side of (5.22), we find 



Kit) 



- 2 



<is 



$i(t, -u z k)e^- Z3)t ^ + $ 3 (t, -uk)e { 

S33 S33 



dk. 



(5.23) 



Using (5.20) and ( 5.23| ) in (5.5a) we find (5.12a). The proof of (5.12b) is similar. 



□ 



Remark 5.3. 1. For the Dirichlet problem, substitution of the expressions (5.11) for 



gi(t) and h\(t) into (5.10) yields a system of quadratically nonlinear integral equations 
for the functions fc)}f. Similarly, for the Neumann problem, substitution of 



the expressions (5.12) for go(t) and ho(t) into (5.10) yields a system of quadratically 
nonlinear integral equations for {&j(t, k)}f . Assuming that these systems have unique 
solutions, S(k) can be determined from (5.8)-(5.9). In fact, we will show in the 
following subsection that these nonlinear systems provide effective characterizations 
of S(k) in the sense that they can be solved recursively to all orders in a perturbative 
scheme. 

2. In the case t hat s ^^jk) has a finite number of simple zeros in S, the repre- 
sentations in (5.11) and (5.12) are valid provided that the integration contour dS is 



replaced with T everywhere, where T is the contour obtained by deforming dS in such 
a way that —V passes above all the zeros of S3s(k) in S, see figure |4j 
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5.4. Effective characterizations. Substituting into the system (5.10) the expan- 
sions 



e 2 $ j2 + 



9o = egoi + e 502 + 



3 = 1,2,3, 
91 = £5n + e 2 5i2 + 



ehni +e h 02 + 



/m 



ehu + e /112 + 



where e > is a small parameter, we find that the terms of O(l) give 3>io 
and $30 = 1. Moreover, the terms of O(e) give $31 = and 



0(e) : *n(t,fc) 
O(e) : *2i(*,fc) 



3 («2-*3)(*-f) 



kgoi + 



mi ]dt>. 



On the other hand, expanding (|5.11|) we find 
3 



2m 
2tH 



ds 



c)S 



2goi 

2a; ^01 
i\/3 



3 

t« ./as 
3w 

7TZ 



(5.24) 
(5.25) 
(5.26) 

$20 = 

(5.27a) 
(5.27b) 

(5.28a) 
(5.28b) 



where S23 = 6^231 + 0(e 2 ) and S13 = 683,31 + 0(e 2 ). The Dirichlet problem can now 
be solved perturbatively as follows. The odd parts of (5.27) yield 

*n_(i, k) = 2u 2 k f e^-^-^hoxdt', $2i-(t, k) = 2k f e^^^ g idt' . 
Jo Jo 

Given 501 and ^01 > we can use these equations to determine $n_ and $21-- We can 
then compute gn and hn from (5.28) and then 3>ii and $21 follow from (5.27). This 
recursive scheme can be continued indefinitely. Indeed, suppose <&ij, $2^ $3j , gu , 
h\j have been determined for all < j < n — 1 for some n > 0. The terms in (5.10) 
of 0(e n ) give 



0(e n ) : $ ln (t 



0(e") : $ 2n (t,k) 



>*) = f 
Jo 



ih 



^0n 



■A J 



dt' + lower order terms, 



3 (z 2 -z 3 )(t-t') 



kgon + "~7=" ) di' + lower order terms, 
v3 / 



0(e n ) : <&3 n (t, k) = lower order terms, (5.29) 
where 'lower order terms' denotes an expression involving known terms of lower order. 



91 



3 
2~Ki 



3oo 
2tH 



k§ 



2gon 



2w/l0n 

iV3 



dk 



Similarly, the terms of 0(e n ) of (5.11) give 

/ + lower order terms, 

./as 

(5.30) 

/ ^-^W^ lower order te rmS . 



3w 



We fi nd $i n _ and <l?2n- from the odd part s of (5.29); then {gi n ,hi n } follow from 
( 5.30[ ) and {$jn}|=i are finally found from (5.29). This shows that for the Dirichlet 
problem {$j}i can be determined to all orders in perturbation theory by solving the 
nonlinear system of Theorem 5.2 recursively. 
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Similarly, for the Neumann problem, the even parts of (5.27) yield 



(5.31) 



while (5.12) yields 



2vr 



h 



01 



as 

2 /o 

2tt 



^ ./as 



□(22-23)* 



'V3 



OS 



7T 



□ (21-23)* 



ns 



(5.32) 



Since g\ and /ii are known, (5.31) can be solved for 3>n+ and $21+) and then 501 can 
be determined from (|5.32 ). Once goi and /iqi have been found, $11 and $21 can be 



computed from (5.27). Extending this procedure to higher orders as in the case of 
the Dirichlet problem, we can determine {^jlf to all orders. 

6. Conclusions 

The 1ST formalism can be used to solve, at least in principle, the Cauchy initial 
value problem for an integrable PDE. However, many physically important problems 
requires solving an IBV problem. In this paper, our goal has been to analyze IBV 
problems for integrable evolution equations with 3x3 Lax pairs. To be concrete, 



we have considered the half- line problem for the system (1.1) whose Lax pair can be 



viewed clS cl 3 x 3 generalization of the Lax pair for NLS. We have shown that: (a) The 
solution can be expressed in terms of the solution of a 3 x 3- matrix RH problem in the 
complex /c-plane formulated in terms of the spectral functions s{k) and S{k). Whereas 
s(k) is defined in terms of the initial data alone, the definition of S(k) requires 
knowledge of all the boundary values, (b) For the linearizable boundary conditions 



(4.1), the function S(k) can be eliminated by using only algebraic manipulation of 



the so-called global relation; the method is therefore as effective as the 1ST on the 
line, (c) In the general case of non-linearizable boundary conditions, the function 
S(k) can be characterized in terms of the known boundary conditions via a system 
of nonlinear integrable equations. This characterization is effective in the sense that 
for small initial data it can be solved in a well-defined perturbative scheme. 

For equations with Lax pairs involving 2x2 matrices, a method announced by 
Fokas in [12j [13] and subsequently developed by several authors has been successful 
in analyzing IBV problems. The results of this paper extend these ideas to equations 
with 3x3 Lax pairs. In fact, even though the transition to 3 x 3 Lax pairs requires a 
number of new steps, the end-result is that the main ingredients of the Fokas method 



can be implemented also for the system (1.1). We expect that our approach will be 



applicable also to other equations with 3x3 Lax pairs. 

Appendix A. Comparison with the case of 2 x 2-matrix Lax pairs 

In this appendix we compare the constructions of this paper with the corresponding 
formalism for 2 x 2-matrix Lax pairs introduced by Fokas, see [14J . According to the 
approach of [14] . the analysis of an integrable PDE on the half- line with a 2 x 2 Lax 
pair of the form 



fix + ifl{k)[(T3, fA = Vifi, 
fit + if2(k)[<J 3 ,l4 = V 2 fi, 



03 







(A.l) 



:-!0 



JONATAN LENELLS 
L3 



proceeds by defining eigenfunctions {[J>j}i as solutions of the Volterra integral equa- 
tions 



lij(x,t,k)=I+ / e-^^W^x' ,t' ,k), j = 1,2,3, 



where 6(x,t,k) = f\(k)x + f 2 (k)t, Wj denotes the closed one-form 



W = e iea3 (V 1 dx + V 2 dt)fi 



(A.2) 



(A.3) 



with ^ replaced by fij, and the contours {7j}f are as in Figure[T] Then, defining sets 
Dj cC,i = l,...,4, by 



L>i = {fc|Im/i > and Im/ 2 > 0}, 
D 2 = {fc|Im/i > and Im/ 2 < 0}, 
D 3 = {fc|Im/i < and Im/ 2 > 0}, 
D 4 = {fc|Im/i < and Im/ 2 < 0}, 



(A.4) 



a RH problem is formulated in terms of the sectionally analytic function M defined 
by 



M = Mr, 



Mi 



M, 




k e D n , n = l, ...,4; 



, [a* 



(A.5) 



M, 



3 1) 



where the spectral functions a,b,A,B are defined b}j^] 
/i 3 (0,0,/c) 



a(fc) 6(fc) 



6(Jfe) a(fc) 



Ati(0,0,fc) 



A(ife) B(jfc) 
Jb® AM 



and = (^ 3 (0,T,jfe))22 = a(fc)A(fc) - b(k)B(k). 

The e igenfunctions {//j}f defined in (2.15) are the 3x3 analogs of the fj,j's defined 
in (A.2). We claim that the eigenfunctions {M n }P used in the main text are the 
3x3 analogs of the M w ' s de fined in (A.5). To see this, we consider the 2x2 analog 
of the integral equation (2.9). Defining \lj,Zj}\ by 

h(k) = -h(k) = -ifi(k), zi(fc) = -z 2 (k) = -if 2 (k), 



we can write the Lax pair (A.l) as 

d(e-^-^ 



W, 



where C = diag(/i, l 2 ) and Z = diag(zi, z 2 ). The -D n 's in (A.4) are given by 

Di = {k\Rel\ > Re/ 2 andRezi > Rez 2 }, 
D 2 = {k\Kel\ > Re/ 2 andRezi < Rez 2 }, 
L>3 = {A;|ReZi < Re/ 2 andRezi > Rez 2 }, 
-D4 = {k\Kel\ < Re/ 2 and Rezi < Rez 2 }. 



We assume, for simplicity, that the symmetries of (A.ll imply that aifjbj(x,t, k)a± = /j,j(x,t,k), 
3 = 1,2,3. 
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(M n )ij(x,t,k) = Sij + / feA^W'WnOr', *',*;)).. 



Consider the eigenfunctions {M n }f defined by the following 2x2 analog of the integral 
equation (|2.9|): 

keD n , i,j = l,2, 
(A.6) 

where the contours 7™-, n = 1, . . . ,4, i, j = 1,2, are defined by (2.10) and W n is given 

by (A. 3) with \i replaced with M n . The following lemma shows that the M n 's defined 
by equation (A.6) coincide with the M n 's defined in (A. 5). 

Lemma A.l. We have M n = M n , n = 1, . . . , 4. 

Proof. Define spectral functions S n (k), n = 1, . . . , 4, by 

S n {k) = M n (0,0,k), keD n . 

Then 

Mn = ^e tx+it S n , n = l,...,4. (A.7) 

Define s(k) and S(k) by s(k) = ^(0,0, k) and S(k) = ni(0,0,k). 

We will give a proof in the case of n = 1; similar arguments apply when n = 2,3, 4. 
The 2x2 matrix ("f 1 )^ '■= is given by 



7 



1 _ ( 73 73 

72 73 



Proceeding as in the proof of Lemma 2.6, we find the following 2x2 analog of the 
matrix factorization problem (2.32): 

s(k;X ) = S 1 (k;X )T^ 1 (k;X ), S{k) = S 1 (k;X )R^ 1 (k;X ), k € D 1 , 

where 



Ri(k;X ) 



Si(k;X 



* * 
* 



Ti(k;X ) 



1 

* 1 



and * denotes an entry yet to be determined. These are eight equations for eight 
unknowns. Solving these equations and letting Xq — > 00, we find 

a{k))' 



Si(k) 

Consequently, by (A.7), 



Hi 



Since the 2x2 analog of equation (2.26) shows that 

be- 2id [ii 2 }i + a[n 2 h = Ma, 



this expression for Mi indeed coincides with the expression for M\ given in (A.5). 



□ 



Remark A. 2. In view of lemma A.l equation (A.6 ) provides an alternative definition 
of the M n 's. In the case of 2 x 2 matrices, since in each domain D n , n = 1, ... ,4, 
there exist two column vectors of the /Uj's which are bounded and analytic, the M n 's 
are conveniently defined by combining the column vectors of the /ij's according to 
(A.5). However, in the case of 3 x 3 matrices, the M n 's involve rather complicated 
combinations of the entries of the /x,-'s, and because of the relatively small domains 
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of boundedness of the //j's, the boundedness properties of these combinations are not 
evident. We were therefore instead led to generalize the alternative definition (A.6). 

We can also verify directly that the 2x2 analog of our definition of the jump 
matrices J m ,n reproduces the jump conditions used in [14J. Indeed, the 2x2 analog 
of definition (2.14) is 

Jm,n = e {S m S n ), 71,171 = 1, . . . , 4. 

For definiteness, we consider the jump across D\ n D 2 . Arguments similar to those 
used in the proof of Lemma A.l show that 

S 2 (k) -- 



^ b(k) 



d(k) 
B(k) 
d(k) 

Thus, the jump condition across D\ n D 2 is given by 



MM a{k) 



M 2 = MiJi 5 2, keD 1 DD 2 where J 1>2 = B (k) 2ie 

\a(k)d(k) e , 

which is the jump condition given in |14j . 



Appendix B. Fredholm integral equations 

We will show that an extension of the standard Fredholm theory can be used to 
establish the required existence and analyticity properties of the solutions {M n }} 2 of 



the integral equations (2.9). This will complete the proof of Proposition 2.2 



Fix n £ {1, . . . , 12} and j £ {1, 2, 3}. Letting Wi(x, t, k) 
write the j'th column of (2.9) as 



(M n )ij(x,t,k), we can 



Wi(x, t, k) = 5. 



f 3 

+ / y~] Kj (x, t; x', t'; k)gwi (x' , t',k), i = 1,2,3, (B.l) 

'«j (=i 



(B.2) 



where the kernel Kj is defined by 

K^t-x'J-^u =e {h-h)(x~x')+^-z 3 ){t-t') 

x (Vi(x' ,t' ,k)dx' + V 2 (x' ,t' ,k)dt')u, i,l = 1,2,3. 

Th ere ar e two reasons why the standard Fredholm theory does not immediately apply 
to (B.l): (a) The equation involves the three different integration contours {7j}f 
instead of being defined on (part of) the real line, (b) The integral kernel Kj is, in 
general, not an L 2 -kernel. Indeed, for k such that ReZj(fc) = Helj(k), the exponential 
factor in (B.2) is bounded, but does not decay as x,x' — > 00. In this case, the kernel 
therefore has decay as x' — > 00, but not necessarily as x — > 00. 

In order to address (a), we first consider the equation on the boundary 



{(x, 0) I < x < 00} U {(0, t) I < t < T}. 



(B.3) 



The idea is to map this boundary to the semi-infinite interval [— T, 00). Define 6{x) 
by 9{x) = 1 for x > and 6{x) = for x < 0. We define a 3 x 3-matrix valued kernel 

JCj(s,s',k) for s,s' G (— T, 00) as follows. Let Kj 



Kf^dx + Kfdt. For s > 0, we 
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let 



JCj(s, s', k) = < 



-Kf(s,0;0,\s'\;k), 



6{s 
0, 



)id x) M; S ',0;fc), 



^ -6{s' -s)K^ x \s,0;s',0;k), 
while, for — T < s < 0, we let 



fCj(s, s , k) 



1 13 

1 1] 

ryV: 
1 13 

ry n . 

1 1] 



0(s' -s)Kf(0,\s\;0,\s'\;k), 
0, 



13 
1 13 
1 13 



-^(O.Hj^O;*), 
Then, introducing the vector valued function v(s, k) by 

'wi(0,\s\,k), -T<s<0 



71 and - T < s' < 0, 
71,2 and < s' < oo, 

72 and — T < s' < 0, 

73 and — T < s' < oo, 

- 7i and — T < s' < cxd, 
= 72,3 and — T < s' < 0, 
= 72 and < s' < oo, 
= 73 and < s' < oo. 



(B.4) 



(B.5) 



Vi(s, k) 



Wi(s, 0,k), < s < oo, 



we can write (B.l) as 

coo 3 



Vi(s,k) 



/OO ' J 
^fCj{s,s\k)uVi{s' 
" T i=i 



, k)ds' , 



-T < s < oo, 



1,2,3. (B.6) 



This is a Fredholm equation of the second kind. It was observed by Caudrey [7] that 
although the kernel JCj is, in general, not of L 2 -type, equations such as (B.6) can be 
analyzed by an extension of the standard Fredholm theory, cf. [18, 30J. Following [7J, 
we find that a sufficient condition on the kernel K,j for equation (B.6) to be solvable 
is that there exists a function m{s' , k) > such that 

\K.j(s, s , k)n\ < m(s , k), i, 1 = 1,2, 3, s, s E (— T, oo), 

and 

/•oo 

m(s',k)ds' = M(k), 



where A4(k) is a bounded function of k. If we consider (B.6) with k G D n , n = 
1, . . . , 12, the boundedness of the exponential factor in (B.2) for k € D n implies 
that such an M{k) exists. We infer that the Fredholm determinant f(k) and minor 
F(s, s' , k) associated with (B.6) are defined and analytic in the interior of the set D n 
and are continuous up to the boundary of D n . Provided that f{k) ^ the solution 
of (B.6) is given by 



Vi(s, k) 



F(s,s',k)ij 



ds . 



-t f(k) 

The zeros of f(k) are the eigenvalues of the Fredholm equation for Vi. Since f(k) is 
analytic this set is discrete. This shows that M n (x,t,k) has the properties stated in 
Proposition 2.2 whenever (x,t) belongs to the boundary (B.3). 

Assuming now that the values of Wi(x,t,k) for (x,t) on the boundary (B.3) are 
known, we can proceed to define M n (x, t, k) also for (x, t) in the set {0 < x < oo, < 
t < T}. Define the vector valued function wo(t, k) by 



w i(t,k) 



6 a + Jo K f ( x > <; °» fc )*/^(o, f, k)df, 75 = 72, 

5 «> 7S=73, 



1,2,3. 
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Then (B.l) can be written as 



poo 3 

(x,t,k) = w (t,k) + / y~] Kj(x, x', t; k)uwi(x', t, k)dx', £ = 1,2, 3, (B.7) 
Jo i=i 



where 



Kj(x, x , t; k)n 



0(x-x')K^(x,t;x',t;k), 7^ = 71,2, 
-9(x' -x)K^\x,t;x',t;k), 7 TJ = 73 . 



For each fixed i £ (0,7"), equation (B.7) is a Fredholm equation to which we may 
apply the approach of [7j. It follows that the solution M n (x,t,k) exists for all (x,t) 



and has the properties stated in Proposition 2.2 
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